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Abstract 

We introduce various optimization schemes for highly accurate calculation of the eigenvalues and the eigen- 
functions of the one-dimensional anharmonic oscillators. We present several methods of analytically fixing 
the nonlinear variational parameter specified by the domain of the trigonometric basis functions. We show 
that the optimized parameter enables us to determine the energy spectrum to an arbitrary accuracy. Also, 
using the harmonic oscillator basis functions, we indicate that the resulting optimal frequency agrees with 
the one obtained by the principle of the minimal sensitivity. 
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1. Introduction 

Eight decades after the discovery of quantum mechanics, the Schrodinger's famous equation still remains 
an interesting subject for various investigations, aiming at extending its applications and at developing 
more efficient analytic and numerical methods for obtaining its energy eigenvalues and stationary states. 
The interest in this subject ranges from various branches of mathematics, physics, and chemistry. This has 
been the driving force behind the development of perturbative and nonperturbative methods for this kind of 
problems. Among them are the factorization method 0,0], semiclassical approximation Q, finite-difference 
technique Q , optimized Rayleigh-Ritz variational scheme [|| , variational matrix solution Q , instanton 
method Q, transfer matrix method Q and many other specific methods. 

One general approach to construct a continuous wave function ip(x) is to represent its values on a set of 
mesh (lattice) points x n which is the starting point of the various mesh methods. Although this approach is 
simple, it is usually very inaccurate which is due to the fact that it only contains the local information of the 
wave function. To overcome this problem, Schwartz proposed a method based on a global construction of an 
approximate wave function which involves only the values ip n at the selected mesh points [Io| . He showed 
that using an optimized mesh spacing the obtained errors are as small as A~ N or even 1/N\ where N is 
the number of the mesh points. This shows the priority of this method over the usual numerical methods 
that yield errors as small as 1/iV, 1/N 2 , etc. The application of a set of orthogonal functions on a finite 
domain for solving a wide class of problems such as function interpolation and the numerical solution of the 



Schrodinger equation has attracted much attention in recent years (see Ref. ll| and references therein) . For 
instance, the variational sine collocation based on the principle of minimal sensitivity (PMS) method can 
be effectively used to obtain the highest precision with a given number of mesh points [l2| . 

The second popular scheme is to expand the wave function in terms of an orthonormal set of the 
eigenfunctions of a Hermitian operator, i.e., the basis-set expansion method. For instance, the trigonometric 
basis functions obeying Dirichlet boundary condition (particle in a box basis) can be effectively used to 
find the spectrum of an unbounded problem. The low lying energy levels are approximately equal to the 
exact ones with high accuracy, if the boundedness parameter is in near vicinity of an optimal value upon 
implementing the Rayleigh-Ritz variational method IH 3] ■ The extension of this method for the periodic 
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boundary condition is also discussed in Ref. [la ]. For these cases, the variational parameter is the width of 
the finite interval (the size of the box). On the other hand, if we expand the wave function in terms of the 
harmonic oscillator eigenfunctions, the frequency of the oscillator would be the variational parameter. 

In this paper, first we briefly outline the Schwartz's method and its error analysis which is based on the 
large distance behavior of the wave function. Then by defining the Schwartz's length, which as we shall 
show has a useful application in diagonalization of the Hamiltonian with the particle in a box basis set, 
the nearly accurate results can be obtained if we use this length as the optimal width of the finite domain. 
By imposing a physically acceptable relation between the potential energy at the optimal length and the 
maximum available basis eigenenergy, we explain the physics behind the Schwartz's length and obtain the 
optimal length without need to the error analysis which is done in the Schwartz's original paper. As it is 
emphasized in Ref. [l^, there is a close connection between the Schwartz's scheme and the Fourier expansion 
where we shall elaborate it in the next section. In this view, we improve our estimation and introduce some 
alternative and more accurate optimal lengths for the particle in a box basis functions. Also we present the 
optimal frequency for the expansion of the wave function in terms of the harmonic oscillator basis functions. 

Another way for finding the optimal length is using the stationarity of the trace of the Hamiltonian Q . 
This method is based on the principle of minimal sensitivity but demands on fixing the values of nonlinear 
parameters before diagonalization of the truncated matrix. In fact, this optimal value extremizes the trace 
and results in highly accurate results. Here we apply this formalism for the trigonometric basis functions 
and find the related optimal length for the anharmonic oscillators. 



2. The Schwartz's method 



Let us consider an analytic reference function u(x) which has simple zeros at the real points 
approximate the wave function iJj(x). We can define an interpolating wave function ip(x) to approximate 

ip(x) as 

u{x) 1 

X " \ Xm ) 



^( X ) = ^Z^ S^i 777Z v W 



ip(x) there. With this definition, the derivatives of ij) at the mesh points — — 

dx 



where ip m — ip(x m ). So at the mesh points x = x n , the interpolating wave function takes the same values as 

or its integrals / tp(x)dx 

only depend on ip m not the derivatives or integrals (see Ref. [10( for details). However, all ij) m would 
contribute to construct the corresponding values of ip which means that the global information about the 
wave function ip is used to find the approximation. 

We can write the above equation as an exact relation by introducing the error term e as 

, / x , / x u(x) 1 

#c)=X>M — U JT-\ +e - W 

X X m U \X m ) 

Note that, we are usually interested for the cases where the wave functions decreases rapidly for large x 
(bound state solutions), so effectively the infinite sum over the mesh points can be truncated to a finite sum. 
In fact, there are two sources of error in our analysis. First one sa due to the analytical approximation 
and the second ex due to the truncation. If we make tA approximately equal to ex by choosing a relation 
between the mesh spacing h and the truncation at n < N, the total error will be reduces considerably. This 
will prevent us from using a too small h when the truncation error dominates or using a too large cutoff 
when the mesh error has the dominant role. 

To compute the errors, consider a bound state wave function which has the following behavior at large 
distances 

ip(x) — e~ ax " , for large x , (3) 

so the truncation error reads 

6T ~ e -a(W (4) 
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For the mesh size error, we need to perform a contour integral in the complex plane. By taking the reference 
function as u(x) = sinOirz/h), the integral can be estimated by the stationary phase method and we obtain 


where q = p/(p — 1) and 



-bh~ 



2(p-l) 



Now the optimal value of h for each TV can be obtained by equating Eqs. @ and 



(5) 
(6) 

(7) 



which results in the exponential decrease of the error by increasing the number of the mesh points 

ewe"™, (8) 

where C = b{a/b) 1 / p . To apply the method, let us consider the following dimensionless time-independent 
one-dimensional Schrodinger equation^ 



-j^+x )4>{x)=E^{x), 



k = 2,4,6,.. 



For this case, the wave function has the asymptotic behavior for large x given by ([3]) with 

k + 2 2 



which results in 
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Now if we define the Schwartz's length L$ = Nhs, we have 



L S (N) 
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(9) 
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where, as we shall show in the next section, it can be used as an accurate candidate for the optimal length 
in the context of the Fourier expansion of the wave function. 

At this point, it is worth to mention the connection between the above collocation method and the 
Fourier expansion scheme. So let us define the generalized sine functions as 



S m (h, x) 



sin [-k(x — mh) /h] 



(13) 



tt(x — mh)/h ' 

where m € Z, uniform grid spacing h and x £R. Now using Eq. ([T]) and u(x) = sm(jrx/h) we obtain 

i>(x) ^^ip m S m {h,x), (14) 



1 Note that for V(x) = j3x k we have E -> /3 k + 2 E. 
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which defines the sine collocation method. It is also possible to write a similar equation in terms of the little 
sine functions Consider an orthonormal set of particle in a box basis functions vanishing at x = ±L 



1 



2L 



(x + L) 



n = l,2, 
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and define 
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where N takes even values. Because of the completeness of the basis functions we have 

N - 

J im 7rr$N(x,y) = S(x - y). 



(16) 



(17) 



By setting h — 2L/N, yk = kh and selecting even values of N, we define the set of (N — 1) little sine 
functions (LSF) as 



s k (h,N,x 



1 fsm[(l + ^) l(x-kh)] cos [(1 + 5^) l{x + kh)] 



2N \ sin [ 2WK ( x ~ kh )] ^[^(x + kh)] 

Therefore, LSF become the standard sine functions when N goes to infinity, i.e., 

/, »r x eva\n(x — kh)/h] r , , 
hm s k {h,N,x) = 1 y ' J = S k (h,x). 

W->oo ir(x — kh)/h 



(18) 



(19) 



The LSF have some common properties with the sine functions, for instance, we can approximate the wave 
function on the interval (— L,L) as 

^(x) = J2^j m s m (h,N,x), (20) 

m 

where can be understood using the definition of Sk(h, N, x) in terms of the completeness relation. Therefore, 
we can rewrite Eq. ([20]) as 



Vm(x). 



In the limit N — > oo this relation becomes 



iP(x) = £ 



ii){x)(p m {x)dx 



-L 



tp m (x) = y^^a m ip m {x), 



(21) 



(22) 



which is the well-known Fourier expansion. So there is a close relation between the sine collocation method 
and the trigonometric basis expansion. 



3. The trigonometric basis expansion 

In this section we study the diagonalization of the Hamiltonian in terms of the particle in a box eigen- 
functions that is basically different from the Schwartz's method. Then by imposing a constraint on the 
potential energy at the optimal length and the maximum available basis eigenenergy, we analytically fix the 
variational parameter. Before going further note that in the Schwartz's scheme N is the number of mesh 
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points whereas in this section N is the number of basis functions. However, in the Schwartz's method the 
reference function vanishes at the N mesh points and the N trigonometric basis functions have at most N 
nodes. Also we have Nh = L. So we expect that there would be a close connection between the Schwartz's 
method and the trigonometric basis expansion which is also explicitly elaborated in the previous section. 
For the potentials which are even functions of x, to avoid large matrices, we can use 

4>in{x) = V7 COS 



\( ^ 


irx 
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T 



/ 1 / 7TL7TX \ 

and 4> m {x) = d — sin y— —J , (23) 



basis functions (to = 1,2,..., N) for even and odd parity solutions, respectively, and write the wave function 
as ip( x ) = Em=i A m 4> m (x) which vanishes at ±L. Now the approximate solutions are the eigenvalues and 
the eigenfunctions of the (N x N) Hamiltonian matrix Hjv where can be written as 



1\ 2 tt 2 . (L 



and 



7T 



H mn = [ to - - J j^5 mn + ( — ) ( Dm+n-i + D m -„ ) , (24) 



H mn T 2 $mn "f - I ) I ^m~n ■&m-\-n I ) (25) 



L 2 

for even and odd states, respectively. Here, S mn is the kronecker's delta and D s is defined as 

2 y, (-iy+ S M fe - 2l - 2 Q 

D s =l I &xx k cos{sx) = \ to [* 2(l+1) (k-2i-l)\ ' (26) 

s = 0. 



k + r 

Now the good strategy is to choose a relation between the basis domain L and the truncation at n < N so 
that the errors due to the basis domain and the truncation approximately suppress each other. This will 
prevent us from choosing a too large domain for small N or a too large cutoff for small L. 

Note that, the expansion the solutions in terms of particle in a box basis functions approximately 
corresponds to confining the potential in an infinite potential wellH i.e., V(x) — x k for |at| < L and V(x) = oo 
elsewhere. Moreover, each energy eigenvalue is a superposition in the form E n — X)m=i ^nm £ m where 
Ylm=i Prim = 1, Pnm > 0, and £ n = n 2 TT 2 / L? . Since for |x| > L this model is not identical with the original 
potential, the basis functions with energies larger than V(L) would not have a useful contribution to the 
sought-after solutions. In this case, we lose accuracy which is due to the small well's width. Also, when 
V(L) 3> £n the solutions would be inaccurate which means that L is too large. Keeping these two points 
in mind, we conclude that the value of the potential at L op should be proportional to the maximum energy 
of the basis functions ejy, namely 

N 2 TT 2 

L k op =a(k) — , (27) 

where a(k) > 1 and it is of order of one. From this equation we can find the optimal value of L as 

i 

L op {N) = (n 2 a{k)^j k+2 N^h, (28) 

which has the same functional form as Ls(N) (TT2"|) . This similarity is due to the relation L — Nh and 
the fact that the maximum number of nodes in this method agrees with the number of mesh points in the 



2 This correspondence is exact as N goes to infinity. 
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Figure 1: The variationally obtained optimal lengths (black diamonds) versus N for V(x) = x k and the predicted optimal 
length curves (120 6 . 



Schwartz's method which both are represented by N . Now, it is only remained to determine a(k) which can 
be found by choosing an acceptable ansatz and comparing with the variationally obtained results (Fig.[T|)[f| 
It is straightforward to check that for the simple harmonic oscillator (k — 2), the most accurate solutions can 
be obtained when the value of the potential at L op is equal to the maximum energy of the basis functions 
ejv, i-e., a(2) = 1. For this case using (|28p we have L op — y/irN — L$- For other values of k, the simplest 

choice is looking for a relation in the form a(k) = rf~z~ with constant r\ to ensure a(2) = 1. As Fig. [T]shows, 
we can properly fit Eq. (1281) to the variational values upon choosing rj = tt/2 which results in 



a° P (*0=(f) 2 . and L op (N) = , ( 29 ) 

Note that, although the coefficients of N^+^ in L$ (fT2^) and L op (|29[) seems to be very different, they are 
approximately equal especially for k < 10. For instance, for k = 4 we have L op — \J n/2 1 / 3 N 1 / 3 ~ 1.579N 1 / 3 
which is nearly equal to L s = 2 1 / 9 n 1 ^N 1 / 3 ~ 1. 5827V 1 / 3 . In Table [Q we have reported the accuracy of 
the energ y ei genvalues for k = {2,4,6,8} which are in complete agreement with the variationally obtained 
solutions [131 ] . 

Since ce(k) for k 1 will not be of 0(1), we expect that the validity of Eq. (|2"9"|) breaks down for large 
k. At this limit we have V(L op ) 3> £n which would result in inaccurate solutions as a consequence of 
too large L op . To check this point, we consider the problem of V(x) = x k where k — » oo. For this case, 
the predicted optimal length is lim/^oo L op (N) ~ ^/7r/2 ~ 1.253. However, the optimal length predicted 
by the Schwartz's formula has the correct limiting value, i.e., lim; c _ ! . 00 L${N) = 1 and gives the following 
proportionality coefficient: 



a s (k) 



\k . 


f-7T\" 


-sm 






Kk) 



(30) 



where as (2) = 1 and as we have desired, it is of the order of unity for all k, i.e., 1 < as < 7r 2 /4. Also, 
Fig. [U shows that the error exponentially decreases with respect to the number of the basis functions. In 
comparison, the calculations show that the usage of Ls results in more accurate energy spectrum than those 
obtained by L op . 



3 From now on we show the relative error of the energy spectrum by e n = 
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Table 1: The relative errors of the energy spectrum of the anharmonic oscillator V(x) = 
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Figure 2: The error versus N for k = 4 and L op (left), and for k = 6 and Lg (right). 



4. Other applications 



The applicability of the introduced optimal lengths is not restricted to the particular form of the poten- 
tials, boundary conditions, or differential equations. Indeed, the optimal length can be effectively used for 
the following issues: 

4-1. Polynomial potentials 

Let us consider a symmetric polynomial potential V(x) = 5Z*L 2 aiX% wnere flfe > and i — 2,4,.... 
Since the optimal length related to each term behaves in a separate manner, we cannot simply use Eq. (|29p 
for this polynomial potential. However, for the large values of L op , the dominant term near the boundaries 
is x . In other words, we can also use this equation for the polynomial potentials when we work with 
a large set of the basis functions. To elaborate this fact, let us study the doubly anharmonic oscillator 
V(x) — a 2 x 2 + a^x A + a$x 6 . The ground state wave function of this potential should not have nodes and it 
should vanish as x — > ±oo. Thus, we can examine the following solution 

V (x) = exp(-^hx 4 + ^b 2 x 2 ^ , b 4 >0. (31) 

It is easy the check that this trial wave function satisfies the corresponding Schrodinger with the eigenvalue 
E Q = —b 2 , b 2 = — 2 a4 °6 ' ^ 4 = a f/ 2 an d a constraint on the potential [l3| 

a 2 = b\ - 36 4 . (32) 

For instance, the potential V(x) — —2x 2 + 2x 4 + x 6 obeys this constraint and have the following ground 
state eigenvalue and eigenfunction 

*o(a0 = exp (~^x 4 - ^x^j , E = 1. (33) 

For this case, using the optimal length related to k = 6, i.e., L op = (7r 2 iV/2) 1 ^ 4 , one can easily obtain the 
ground state energy E with 28 and 36 significant digits accuracy for N — 30 and N — 40 basis functions, 
respectively. 

4-2. Periodic boundary condition 

We have deduced Eq. ([2T))) from diagonalization scheme for the Hamiltonian with the basis functions 
which are vanishing at the boundaries, but we can also use it for the case of the periodic boundary condition 
under some circumstances. Note that, for this case, we observe an inflection point instead of a minimum 
value in the graph of the energy versus the domain of basis functions (ljj . However, since for the large 
optimal lengths the low- lying wave functions are almost zero at the boundaries, Eq. (|29[) is still valid for 
case of the periodic boundary condition when L op is large enough. 

4-3. Multidimensional problems 

In a d-dimensional space which the Hamiltonian is invariant under parity transformation, the basis 
functions are in the form 

c^-m.^ = nv? 008 ^-^)^!' ( 34 ) 

<U~m,(2) = n/fsm(^V (35) 

i— 1 ^ ' 

where m t = 1, 2, . . . , N. Now the Hamiltonian is expressed as a N d x N d matrix. Therefore, even for the 
small values of N, the resulting matrix is too large and the ordinary variational scheme cannot be applied 
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efficiently. In general, we need to find d optimal lengths for the accurate calculation of the eigenvalues 
and the eigenfunctions. For instance, for the separable potential V(x) — Yli=iV( x i)i we can use the 
introduced optimal length (12T)1) which needs to be properly chosen for each direction to get the desired 
accuracy. Moreover, if the Hamiltonian has the rotational symmetry, i.e., V(x) = we only need one 

optimal length to diagonalize the N x N Hamiltonian. 

4-4- Wheeler- De Witt equation 

One approach to quantize gravity is based on the well-known Whcclcr-DcWitt equation (WDW), i.e., 
= where ^ is the wave function of the universe [lfjj]. After freezing out many degrees of freedom, 
WDW equation is expressed in the minisuperspace and ^ is a function of just a few variables such as the 
scale factor, the scalar field, etc. This equation is a hyperbolic partial differential equation and can be cast 
in the form of the Schrodinger-like equation in the vicinity of the boundaries. To solve the WDW equation, 
the proper adjusting the width of the domain of the basis functions with the number of the basis functions 
is also an important issue which results in accurate solutions 17- 19|. 



5. Alternative proposals 

To improve the estimation for the optimal length, we need to modify Eq. (|29[) in such a way that it takes 
a finite value as k goes to infinity. For instance, consider the following proposal 

k ~ 2 /„4(&— 1) \ fc(fc + 2) 

^(*0=GEr, -d 4 1 p ) (iV)=(^U 1 ) N*. (36) 

For k = 2 and k = 4 this proposal coincides with Eq. (f2l?l) and for large k it agrees with Eq. (|3H)) . i.e., 

Jim <*#(*)= lim a s (k) = ^-. (37) 

Since a p is nearly equal to as for all k (see Fig. [3J, the accuracy of this proposal is of the order of as- 

Up to now, the most accurate proposal for the optimal length is Ls- However, the calculations show that 
to find more accurate results, we need to slightly increase Ls for all k and N. Therefore, since we always 
have Lop < Ls, the accuracy of Eq. (|3"6"j) is slightly lesser than Eq. (|2l?)) which can be also confirmed by 
explicit calculations. To find an optimal length that is more accurate than Ls, note that the factor it/2 in 

k / 7T \ 

Eq. (|36[) is the first term in the asymptotic expansion of the term — sin ( — ) in Eq. (|30p . i.e., 

2 V k J 



dn (D=f-w + °^ (38) 



So a good idea is to add the second factor ^|p- to it J 2 in Eq. ((36)) and write the optimal length as 



fc-2 fc— 2 

fe/2 / _2 \ 7T/2 



and 



1 2(k — 2) 2(fc — 2) 

/ 4(k-l)\ k(k+2) / 2 \ k(k+2) ( n 2 \ fc( fc+ 2) 

L -^=(W) ( 1+ 6^) iV ^ = ( 1+ 6^) ^ )(7V) ' (40) 

which satisfies a op (2) = 1 and lim^oo a„ p !(k) = 7r 2 /4. It is straightforward to check that L op (a op ,) 
is always slightly greater than Ls (as) (see Fig. [3]). Indeed, calculations show that between the several 
mentioned candidates for the optimal length, L op gives the most accurate energy spectrum. As a simple 
application, in Table [5J we have reported the ground state energy of the anharmonic oscillators using only 

(2) 

one basis function (N = 1) and the optimal length L op . For these cases, the matrix of the Hamiltonian has 
only one element and the relative error is less than 20%. 
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Table 2: The ground state energies and their relative errors using one basis function (N = 1) and L p 




Figure 3: The proportionality coefficient a(k). 



6. Stationarity of the trace of the Hamiltonian 

In this section we apply the optimized Rayleigh-Ritz method proposed by Okopinska 0, Q for the 
trigonometric basis functions. This method is originally used for the harmonic oscillator eigenfunctions 
and it is based on fixing the values of nonlinear parameters before diagonalization of the truncated matrix. 
Before diagonalization, the only physical quantity that can be determined is the trace of the Hamiltonian 



N 

Tr N H = ^(n\H\n) 

n=l 



(41) 



which represents the iVth-order approximation to the sum of energies of the N lowest bound states. Now 
the strategy is to choose the the values of nonlinear parameters so as to make Tr^H stationary. Since the 
only nonlinear parameter for the trigonometric basis functions is L we have 



dL 



Tr N H = 0. 



Using Eq. (|2~4"|) we can find Tr^iJ as 



Tt n H 
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where D = ir k /(k + 1) and 
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Figure 4: The optimal lengths versus N for k = 2 (left) and k = 4 (right). 



where £(s) and £(s, a) are the Riemann zeta function and the generalized Riemann zeta function, respectively. 
After some algebra, Eq. (|4"2"j) results in the following optimal length and «t 
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a T (fc) 
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3k 
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For fc = 2 we have L^(N) ~ \/ttN which agrees well with L$ and L p- However, the accuracy of this method 
reduces for k > 2. In comparison with L op , we should mention that for k < 10, L op is more accurate as it 



is apparent from Fig. 2] and l.h.s of Fig. [5o But for k > 10, £t is more closer to L$ (r.h.s of Fig. [S]) and 
similar to tends to one as k goes to infinity. The accuracy of the energy spectrum of the anharmonic 
oscillators for k = {2,4,6,8} and N = 10 is reported in Table [3] As the table shows, the optimal length 
obtained using Ls or L op formulas gives more accurate results in comparison with Lt formula for a given 
N. In Fig. [SI we depicted the coefficient of proportionality a{k) for the proposed schemes. Since a op grows 

exponentially, it cannot be used efficiently for large k. On the other hand, as, 0$ , and aop* go to 7r 2 /4 and 
«t goes to 2/3 at this limit. 



7. The harmonic oscillator basis expansion 

We can also find the energy spectrum of the anharmonic oscillators using the basis of the harmonic 
oscillator eigenfunctions 
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4 Note that for k = 2, we have L op = Ls, and for k = 4, L op nearly coincides with Ls- 
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Figure 5: The optimal lengths versus N for k = 6 (left) and A; = 16 (right). 
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Figure 6: The proportionality coefficient versus k. 
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Table 3: The relative errors of various optimized schemes for N = 10 basis functions. 
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where H n (flx) are Hermite polynomials and the frequency Q is the variational parameter. This frequency 
can be fixed by the principle of minimal sensitivity, requiring the dependence on the variational parameter 
be as weak as possible [20j. For the Hamiltonian 



H = l 
2 



d 2 , 2 2 



Aa; 4 , (48) 
the application of the PMS to the sum of N even and N odd basis functions gives [H[ 



^ MS - w 2 ttpM S =8A(V+^. (49) 

It is also possible to find a similar relation for the optimal frequency fi p using the prescription presented 
in Sec. [3] Here we have two potentials: one physical ^u! 2 x 2 + Xx 4 and one unphysical \Sl 2 x 2 which we use 
the eigenfunctions of the latter to approximate the solution of the former. The intersection points of these 

potentials are given by x- m t — n2 ~j* ^-. So the potentials take the following value at these points: 

V(x illt ) = ^(n 2 -cj 2 ). (50) 

On the other hand, for N even and N odd basis functions, the maximal quantum number is 2N — 1 and 
therefore the maximal energy reads 

fimax =n\2N-^j. (51) 

Now since the basis functions with the energy larger than V(xint) would have no useful contribution to the 
sought-after solutions, by equating Eq. (1501) and Eq. (|5ip we obtain the relation for the optimal frequency 
as 

Q 3 op -uj 2 Q op = 8x(n-^\ , (52) 

which agrees well with Eq. (j4"9"]l especially for large N. At this limit, we have Q op ~ fipsM ~ 2(A7V) 1 / 3 . 
Also, the N dependence of the optimal frequency coincides with the N dependence of the optimal length 
in the trigonometric expansion for k = 4, i.e., Q op ~ N 1 ^ 3 ~ L op . Note that, the accuracy of f2 p is of the 
order of SIpsm even for small N. In Fig. [JJ we have depicted in terms of the total number of the basis 

functions 2N for uj 2 = {0, 12A 2 / 3 , 20A 2 / 3 }. 



8. Conclusions 

In this paper, we presented several optimal lengths for accurate calculation of the eigenvalues and eigen- 
functions of the anharmonic oscillators using trigonometric basis functions obeying Dirichlet boundary con- 
dition and the harmonic oscillator eigenfunctions. We indicated that the value of the potential at the optimal 
length should be proportional to the maximum energy of the used basis functions. Since both V(L op ) > ejv 
and V(L op ) <C £jv are two sources of error, we demanded that the proportionality coefficient a{k) to be 
of the order of one. For the trigonometric basis set, we suggested some ansatz for a{k) with the condition 

(2) 

a(2) = 1 and found that L op gives the most accurate results. By defining L$ we showed that it can be 
also used as an accurate optimal length. Indeed, it is shown that the optimal number of the mesh points 
(for fixed h) in the Schwartz's scheme, where the reference function vanishes there nearly coincides with 
the optimal number of the basis functions (for fixed L) which at most have the same number of nodes. 
Moreover, as is of the order of unity for all k. An alternative proposal is using the trace of the Hamiltonian 
as an only physical quantity before diagonalization and to make it stationary at the optimal value L-r- We 



13 




Figure 7: The optimal frequency parameter versus 27V obtained via the principle of minimal sensitivity for ui 2 = 
{0, 12A 2 / 3 ,20A 2 / 3 }. 



showed that this optimal length has the correct asymptotic value for large k and can be used as a good 
approximation for k > 10. For all proposals except L op the proportionality coefficient a(k) remains of order 
of one for all k and they have the correct asymptotic value, i.e., limfe_>. c>0 L(N) — 1. We indicated that these 
proposals can be also used for the multidimensional problems, and for large N, for the periodic boundary 
condition and the polynomial potentials. The following schematic diagram shows the relative efficiency of 
the proposed optimal lengths 



4p } . (53) 



where the right arrow indicates the direction of the increasing of the accuracy. For the harmonic oscillator 
basis functions, we showed that the resulting optimal frequency agrees with the one obtained using the 
principle of the minimal sensitivity. 
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